Analyses of multiplicity distributions by means of the 
Modified Negative Binomial Distribution and its 

KNO scaling function 

Takeshi Osada, Noriaki Nakajima^, Minoru Biyajima^ 
and Naomichi Suzuki^ 
Department of Physics, Tohoku University, Sendai 980-77, Japan 
■^RCNP, Osaka University, Ibaragi 567, Japan 
^Department of Physics, Shinshu University, Matsumoto 390, Japan 
^Matsusho Gakuen Junior College, Matsumoto 390-12, Japan 

February 1, 2008 

Abstract 

We analyze various data of multiplicity distributions by means of the Modified 
Negative Binomial Distribution (MNBD) and its KNO scaling function, since this 
MNBD explains the oscillating behavior of the cumulant moment observed in e~^e~ 
annihilations, h-h collisions and e-p collisions. In the present analyses, we find that 
the MNBD (discrete distributions) describes the data of charged particles in e~^e~ 
annihilations much better than the Negative Binomial Distribution (NBD). 

To investigate stochastic property of the MNBD, we derive the KNO scaling 
function from the discrete distribution by using a straightforward method and the 
Poisson transform. It is a new KNO function expressed by the Laguerre polyno- 
mials. In analyses of the data by using the KNO scaling function, we find that the 
MNBD describes the data better than the gamma function. 

Thus, it can be said that the MNBD is one of useful formulas as well as NBD. 



1 Introduction 

Recently it has been found that cumulant moments of the multiplicity distributions 
in both e~^e~ annihilations and hadronic collisions show prominent oscillatory behavior 
when plotted as a function of their order g Jll] . In Ref . [jl| this behavior is attributed to the 
Q CD-type of branching processes apparently taking place in those reactions. However, 
in Refs. p|, ^ we have shown that the same behavior of the moments emerges essentially 
from the modified negative binomial distribution (MNBD) (which actually describes the 
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data of negatively charged particles much better than the negative binomial distribution 
(NBD) 0, 1^, §). This distribution can be derived from the pure birth (PB) process with 
an initial condition given by the binomial distribution [0]^. 

In this paper, first of all, we analyze various experimental data of charged particles 
20|] including e-p collisions ||2l]] by the MNBD and the NBD, to elucidate why the 



MNBD describes the data better than the latter in e^e~ annihilations 0, ^. 

Second, we derive the KNO scaling function of the MNBD both by the straightfor- 
ward method (i.e., proceeding to the limit of large multiplicities n and large average 
multiplicities (n) while keeping the scaling variable z = n/{n) finite and fixed) and 
by the Poisson transform. Using this KNO scaling function we analyze the observed 
multiplicity distributions in e+e~ annihilations, h-h collisions and e-p collisions. 

Finally the concluding remarks are given. The generalized MNBD is discussed in 
Appendix. 



2 Analyses of data by discrete distributions 

It is known that the MNBD is obtained from the following equation governing the PB 
stochastic process^, |^, |^, |^, 

dP{n, t) /dt = -XnP{n, t) + X{n - l)P{n - 1, t) (1) 

with an initial condition 

P(n, t = 0)= - a)^"", (0 < a < 1). (2) 

Here A is a birth rate of particles. The parameter t describes the evolution of the 
branching processes. The generating function (GF) of the distribution at t = is given 
by 
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n(M; t = 0) = ^ P(n, t = 0) X 



n=0 
oo 



N-n^n 



n=0 



= [au+l-a]^, (3) 
The GF for the MNBD at t = f^dt = T is given as Refs.||, |, 0, 



oo 

n{u;T) = Y.P{n,T) x u"" 



n=0 



^ It is interesting to mention here that it comes also from the concept of purely bosonic sources as 
presented recently in Ref.H. 
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(4) 



where 



-1 — ri(M — 1)- 
-1 - r2{u - 1)- 

ri = -«(! +p) +p, 
r2=p (> 0), 
p = exp(AT) — 1, 

and N (an integer) corresponds to the number of possible excited hadrons at the initial 
stage. The term AT is evaluated by AT = ln(l + r2). The MNBD is given by the GF, 
Eq.(4), as 



P(0) 



Pin) 



n(o) 

1 d'^Uiu 



1 + r2 
1 



n\ du"- 



nl 



N 



r2 



1 + r2 



r2 - ri 



;i + r2)^ 



(5) 



ri 



r2 



Co 



-a 



1 

N 



{n) 



1 

2' 



1 + 1^ (n) - -. 



Equation (||) is applied to various experimental data[P]-[pl[]. 

For the sake of comparison with results obtained by Eq.(^), we also use the negative 
binomial distribution (NBD) given by 
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T{n + k) ( k 



T{n + l)T{k) \{n) J 
where k > 0. The GF of the NBD is given by 

(") 



-(n+fc) 
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k 



(6) 



(7) 



The GF, Eq.(4), of the MNBD reduces to that of the NBD, if n = 0. 

Both the MNBD and the NBD are applied to analyses of the experimental data in 
e~^e~ annihilations|0-|TB|, h-h collisions ||17||-[P(]|] and e-p collisions ||21||. Since there are 
maximum values of multiplicity observed, n^axi we introduce a possible bound for 
and truncate the multiplicity distribution P{n) at rtmax and renormalize it as follows. 



C X ^ P{n) = 1 

n=0 
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In analyses in Ref.[^ 0, authors used the following treatments; 1) In the first case, the 
data (n) is used, ri and are free. Only a constraint, r2 = ri + {n)/N, is used. 2) In 
the second case, (n), ri and are free. r2 = ri + {n)/N is used. 

On the other hand we use C2 and (n) as inputs. Only is free. Our results by 
means of Eqs.(^) and are given in Table I and some typical results of fitting are also 
given in Figs. 1(a) and (b). It is found that in e+e~ annihilations, the minimum values 
of x^'s obtained by fitting Eq.(^) to the data are much smaller than those of the NBD. 
This fact corresponds to our previous work in which we have found that the Hg = Kg/ Fg 
moments obtained from the MNBD are much better than those of the NBD in e^e~ 
annihilation[2,3]: However in h-h collisions they are almost equivalent. 

Because the data for e-p collisions by HI Collaborations have recently reported, 
we calculate Hg moments by the MNBD and the NBD. Two results of comparisons 
are shown in Figs. 2(a) and (b) for energies -y/i = 115-150 GeV and y/s = 150-185 GeV, 
respectively. It is found that the results of the MNBD and the NBD are almost equivalent 
as observed in h-h collisions 0. 



3 The KNO scaling function 
3.1 The straightforward method 

In order to know stochastic property of the MNBD, we consider the KNO scaling func- 
tion of Eq.(^. Traditionally the KNO scaling function is derived from the multiplicity 
distribution P{n) multiplied by the corresponding mean multiplicity (n) by going to 
the large multiplicity n and large mean multiplicity (n) limit while keeping their ratio, 
z = lim„^(„)^oo "^Z {n) fixed. In our case, starting from Eq.(|^) we arrive at the following 
function 

"^{z) = lim {n)P{n) 

n,(n)—> 00 

= 0'.-*-|.c.^(-^)'(M)'..-, ,, 



The parameters and r'2 in Eq.(^ are given by 

^ 2 





1 


1 - 




iV 






1 + 




iV, 



\{c2-l + l-){n), (9) 



which are slightly different from ri, r2 in Eq.(^, because {n) ^ 1. It should be noticed 
that the normalization of Eq. (H) differs from the unity, 

f = l- M . (10) 
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where the second term corresponds to the term (n)P(O) in Eq.(|^). 



3.2 The Poisson transform 

The KNO scahng function t) is related to the muhiphcity distribution function 
P{n,t) by the Poisson transform 



, , inverse Poisson trans. ^ , , 

Poisson trans. 



In this approach we obtain that 



P{n, t) = I -r_^e— VP (j^, t] , (11) 

JO nl \{n)/a J {n)/a 



{n)/a J 271 a J-oo V a 



where Il{u,t) is the generating function of P{n,t). These equations hold also at t = 
(the stationary function). 

Using the generating function Eq.(^), the KNO scaling function "^{z) is given by the 
following inverse Laplace transform 

1 p^.^(ri)"f:„c,(!i^My7, + M)-',,. (13) 

Then we arrive at the KNO scaling function for the MNBD, 



A J V2J ~{ r(j) \ r[ J \ r. 



'2 



where L^"\x) is the associated Laguerre's polynomial. In Eq.(0) the first term corre- 
sponds to the constant term in Eq.([T3|) (, or (n)P(O)). In the numerical calculations 
the first term is very small because {r'i/r'2)^ where (A^ > 1). Therefore Eqs.(7) and 
(13) are almost equivalent in numerical analyses. The generalized MNBD is discussed in 
Appendix. 
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The Cq moment of the MNBD is given by 



roo 



Jo 



j=i ^ vj; \ '1 / \ ' 2 




We also analyze the data by the gamma function which is the KNO scaling function 
obtained from Eq.(^); 

^^^^{z) = ^f'-'^'-'. (16) 
3.3 Analysis of the experimental data 

We investigate the applicability of the MNBD as presented in its KNO form (i.e., 
using Eq.(|I^)) to the description of the observed multiplicity distributions in e+e" 



annihilations P| - 1|1 6|| , h-h collisions p^, pO|] and e-p collisions ||2T]]. In actual analysis the 
values of the r'^ and r'2 are determined by using the experimental data of C2 and (n) 
as menstioned before. Table II shows obtained parameters in Eq.([T^) and minimum 
values of x^'s of fitting to the experimental data. See Figs. 3 (a) and (b). We show here 
some typical results of fittings. The results using Eq.(|16D are also shown there. As is seen 
from Table II, the minimum values of x^'s for the MNBD fitting are a little smaller than 
those for the NBD fitting in low energies e^e~ annihilations and h-h collisions (below 
\/s ~ 50 GeV). For e-p collisions the minimum values of x^'s are also smaller than those 
of the NBD up to energy v^=220 GeV. 

On the other hand, in high-energy (LEP and SPS energies for the e+e~ annihilations 
and h-h collisions, respectively), values of the minimum x^'s by the MNBD fitting are 
almost equivalent to those of the NBD fitting. 



4 Concluding remarks 

Through analyses of the data by discrete distributions, we find that the MNBD describes 
the data in e^e~ annihilation much better than the NBD. See Table I. We find that 
iV = 7 ~ 13 except for the data by HRS and TASSO Collabs. On the other hand, 
in |] A^ = 7 ~ 8 for negatively charged particles is obtained, provided that their 
treatments for parameters are used. In other words, we have to pay our attention to the 
determination of parameters contained in the formulas. 
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Second the Hq moments of e-p collisions by HI Collab. are analyzed by two distribu- 
tions. The results are almost the same as h-h collisions. The 's in Table I are slightly 
better than those in Refs.|^ when our treatment for parameters is used. 

Third, the KNO scaling function of the MNBD is obtained. As far as our knowledge 
for the KNO scaling functions [p^ ] is correct, we have not known the KNO scaling function 
expressed by the Laguerre polynomials. It is applied to the analyses of the observed 
multiplicity distribution in e+e~ annihilations, h-h collisions, and e-p collisions. The 
data are also analyzed by the gamma distribution. As is seen from Table II, the MNBD 
describes the experimental data better (for e^e~ annihilations) than or as well as (for 
h-h and e-p collisions ) the NED. 

In conclusion, it can be stressed that the MNBD is useful as well as the NED. See 



also Ref. 25 



What we have found in the present analyses suggests us the following: In e+e~ 
annihilation, the stochastic pure birth process with the binomial distribution at t=0 is 
useful, becasuse the finite number (corresponding to = 7 ~ 13) of the excited hadrons 
(or the pair creation of quarks) is probably expected. In other words, the binomial 
distribution as the initial condition is more realistic than other condition {5n,k) in the 
stochastic approach. 

To know stochastic property of the MNBD, we have considered the generalized 
MNBD. In its concrete application, we have known that the discrete distribution of 
Eq.(A5) can not explain the oscillating behaviors of the cumulant moment observed 
in e+e~ annihilations and in hadronic collisions much better than the MNBD, Eq.(^). 
For example, x^'s values for discrete analysis of data at a/s=900 GeV are as follows: 
(Ar,fc,x^) = (l, 2.22, 65.2), (2, 1.33, 63.0) and (3, 0.06, 60.8). These results suggest us 
that it is difficult to determine the best combination of parameters, because of three pa- 
rameters. It seems to be necessary more carefulness and skillfulness than the use of the 
MNBD. Thus the generalized MNBD is given in Appendix. However, we are expecting 
that Eq.(A5) will become useful in analyses of data of some reactions at higher energies, 
since it has both stochastic characteristics of the MNBD and the NBD. 
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Appendix: 

In order to know the stochastic structure of the MNBD in detail, we discuss the following 
point: The solution obtained from the branching equation of the pure birth process with 
the immigration under the initial condition of the binomial distribution]^, |^ is one of 



7 



the extensions of both the MNBD and the NBD. In this case the stochastic equation 
(0) should be changed as follows: 



dP{n, t) /dt = -\nP{n, t) + \{n - l)P{n - 1, t) - AoP(n, t) + \oP{n - 1, t), (Al) 
where Aq is an immigration rate. Its generating function is given as 



n(«) = ^ P( 



n]u 



[l-f,(u-l)f[l-f2{u-l)] 



-k-N 
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and k = Aq/A. The MNBD is obtained by neglecting the power k in Il{u). The gener- 
ating function of the NBD is given by neglecting the power in Eq.(A2). The physical 
meaning of the immigration term k may be interpreted as a possible contribution from 
constituent quarks and gluons. Using Eq.(0), we have directly the KNO scaling function 
for Eq.(A2) 
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(A5) 



Here the term 1 / (n) in the parenthesis in the squared root of Eq.(A3) is neglected because 
(n) — >■ oo. It also noticed that all factors (n)'s in Eq.(A5) are canceled out. This function 
becomes the KNO scaling function of the MNBD Eq.(p!4D when /c — 0, and reduces to 
the gamma distribution Eq.(|l6|) if A^ = 0. 
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Table captions 
Table I 

The minimum values of x^'s for fitting to the experimental data using the MNBD and 
the NBD, respectively. The data of (n) and C2 are used. We also show giving the 
minimum values of x^'s in Eq.(6). In analysis of the data by UA5 CoUab (y^=900 GeV 
), we find a minimum at A^=3 giving Pin) < at small integer n. We discard this 
and adopt = 4 giving Pin) > for all n. In analysis of the data by DELPHI Collab. 
(n) = 21.26 and C2 = 1.0908 are calculated from the data. 

The symbole * denotes the minimum value of x^'s between 2 columns, x^'s in Ref.|0] 
are obtained by means of the statistical errors only. 



Table II 

The same as Table I but for fitting of the KNO scaling function obtained from the MNBD 
and the NBD, respectively. The data of (n) and C2 are used. In analysis of the data 
by OPAL (v^=133 GeV ) and UA5 {^/s =900 GeV) Collabs, we find minimum x^'s at 
A^ = 8 and A^ = 3 giving P{n) < in small z, respectively. We discard these and adopt 
A^ = 10 and A^ = 4 giving P{n) > for all n for the data by OPAL and UA5 Collabs, 
respectively. 

The symbol * denotes the minimum value of x^'s between 2 columns. 



Figure captions 

[Fig. 1] Results obtained by discrete distributions. 

(a) The result of fitting by Eqs.(5) and (7) to the data observed by TASSO Collab [[TO 
(Vs=34.8 GeV.) 

(b) The same as (a) but for the data observed by OPAL Collab |jl6|]. (y^=91.2 GeV.) 

[Fig. 2(a) and (b)] Analyses of Hg moment in e-p collisions by HI Collaboration for 
energies (a) ^5=1 15-150 GeV and (b) v^= 150-185 GeV. 

[Fig. 3(a) and (b)] Results obtained by the KNO scaling functions. 
The same as Figs. 1(a) and (b) but for fittings Eqs.(P^ and ([TB|). 
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1 1 1 1 


1 1 1 


HI 115-150 GeV 


- 






1 <^ n* ^ ^ 
1 ^ 1 1 ^ J 








© Exp 


- 


1 




NBD 








MNBD 


- 












• 




- 








- 


• 




- 




, , , , 


, , , 


with truncation 

1 1 1 1 1 



5 10 15 




